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Self-trapping  of  beams  in  photorefractive  semiconductors  offers  several  attractive 
features  when  compared  with  more  conventional  materials  such  as  strontium  barium  niobate.  For 
example,  operation  at  optical  telecommunications  wavelengths  and  fast  response  times  combined 
with  very  low  operation  power  levels  are  a  few  of  the  advantages  offered.  These  properties  can 
lead  to  many  exciting  applications,  such  as,  reconfigurable  switches,  interconnects,  and  self- 
induced  waveguide  devices.  Unfortunately,  photorefractive  semiconductors  tend  to  have  tiny 
electro-optic  coefficients,  which  imply  that  a  "conventional"  screening  soliton  in  such  materials 
would  require  applied  fields  that  are  too  large  for  most  applications.  But  Fe-doped 


photorefractive  InP  crystals  displays 
photorefractive  effects  that  stem  from  excitation 
of  both  electrons  and  holes,  and,  when  the 
photo-excitation  rate  of  holes  is  comparable  with 
the  thermal  excitation  rate  of  electrons,  the 
induced  space  charge  field  is  resonantly- 
enhanced  by  more  than  an  order  of  magnitude. 
That  is,  the  internal  photo-induced  field  exceeds 
50  kV/cm  even  when  the  applied  (bias)  field  is  5 
kV/cm  or  less.  This  enhancement  appears  only  at 
the  vicinity  of  a  particular  intensity  of  the  beam. 
Furthermore,  if  a  given  polarity  of  the  applied 
field  leads  to  self-focusing  at  intensities  below 
the  resonant  intensity,  the  same  polarity  leads  to 
beam  breakup  and  defocusing  when  the  intensity 
exceeds  the  resonant  intensity.  These  features 
are  truly  unique  to  photorefractive  solitons  in 
InP  and  cannot  be  explained  by  the  established 
theory  of  photorefractive  screening  solitons.  In 


Fig.  1.  One-dimensional  interconnects  in  InP 
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fact,  all  the  theories  describing  all  soliton  effects  in  photorefractive  do  not  show  any  resonance  of 
any  sort. 

A  photorefractive  soliton''^  is  created  when  a  photo-induced  index  change  exactly 
compensates  for  the  diffraction  of  the  beam.  In  this  sense  the  beam  is  able  to  create  its  own 
waveguide  and  offers  potential  for  applications  in  the  field  of  all-optical  switching  and  beam 
steering.  These  effects  have  been  extensively  studied  in  ferroelectric  oxide  and  sillenite  oxide 
crystals  for  visible  wavelengths.  For  the  near  infra-red  wavelengths  used  in  telecommunications, 
InP:Fe  crystals  have  already  demonstrated  interesting  photorefractive  properties and  self¬ 
trapping  of  a  one-dimensional  beam  had  been  reported  by  our  research  group  (fig.l).‘‘‘  In  the 
research  supported  by  this  grant  we  report  on  the  first  observation  of  trapping  of  a  two- 
dimensional  Gaussian  beam  at  1.3  pm,  guiding  of  a  second  beam  at  the  telecommunications 
wavelength  1.55pm  and  on  the  results  from  collisions  between  two  spatial  solitons. 

EXPERIMENTAL  RESULTS 

For  the  experiment  on  2-D  waveguide  formation  the  beam  from  a  Nd:  Yag  laser  at  1 .3pm 
is  collimated  and  focused  with  a  5cm  focal  length  lens  on  the  entrance  face  of  an  InP:Fe  crystal 
whose  temperature  is  stabilized  at  297K.  The  electric  field,  Eo,  is  applied  along  <1 10>,  while 
the  beam  propagates  along  <1  1  0>,  and  is  polarized  either  horizontally  or  vertically  at  45°  from 
<110>. 

Figure  2  shows  an  example  of  2-D 
trapping  using  this  configuration.  The 
beam  size  at  the  entrance  face  of  the 
crystal  has  a  diameter  of  about  55pm 
and  diverges  to  a  170pm  diameter  at  the 
exit  face  when  no  field  is  applied  to  the 
crystal.  The  InP:Fe  crystal  length  is  1cm 
in  the  direction  of  the  propagation. 

When  a  12  kV/cm  field  is  applied  to  the 
crystal,  the  beam  is  trapped  and  the 
beam  diameter  at  the  exit  face  is 
elliptical  and  reduced  to  about  55pm  by 
65  pm.  Our  experiments  show  that  for 
exactly  the  same  value  of  all  parameters, 
the  beam  is  trapped  to  the  same  diameter  for  either  a  1cm  or  a  0.5cm  crystal  length,  supporting 
the  conclusion  that  the  trapped  diameter  is  the  same  throughout  the  crystal.  To  show  that  an 
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Fig.  2.  Two-dimensional  interconnect  formation  in  InP 
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efficient  waveguide  is  formed  in  the  crystal,  we  have  also  propagated  a  1.55pm  (=  ImW/cm^) 
laser  beam  collinear  to  the  1.3pm(lW/cm^)  trapped  laser  beam.  To  observe  only  the  1.55pm 
beam,  a  color  filter  is  placed  after  the  crystal.  When  both  beams  are  present  and  vertically 
polarized,  the  1.55pm  beam  is  observed  to  be  guided  along  the  same  direction  as  the  1.3pm 
beam.  However,  when  the  polarization  of  the  1.3pm  beam  is  changed  to  horizontal  and  the 
1.55pm  beam  kept  vertical,  the  1.55pm  beam  is  guided  on  the  side  of  the  1.3pm  trapped  beam. 
While  in  both  cases  the  1.55pm  beam  is  guided,  we  can  see  that  when  the  two  beams  have 
orthogonal  polarizations,  the  1.55  pm  beam  is  more  efficiently  guided  and  has  a  near  circular 
profile.  The  background  noise  is  at  least  an  order  of  magnitude  below  the  trapped  beam. 

In  order  to  measure  the  index  change  responsible  for  this  large  focusing  effect,  an 
interferometer  is  added  to  the  apparatus  (Fig.  1).  Specifically,  a  1.3pm  wavelength  probe  beam 
is  launched  through  the  crystal  to  monitor  the  refractive  index.  The  probe  beam  is  an  expanded 
beam  polarized  perpendicularly  to  the  focused  beam.  It  is  also  lower  in  intensity  than  the 
focused  beam,  so  that  it  does  not  perturb  the  photorefractive  space-charge  field.  With  this 
orientation  an  index  change  with  the  same  magnitude,  but  with  opposite  signs,  is  seen  by  the  two 
beams.  A  rotatable  polarizer,  positioned  after  the  crystal,  selectively  isolated  either  beam  for 
imaging  onto  the  camera.  When  the  probe  beam  is  imaged  on  the  camera,  a  reference  beam 
which  is  identically  polarized  is  also  present.  In  this  way  an  interference  pattern  depicting  the 
index  change  in  the  InP:Fe  crystal  is  obtained. 

To  evaluate  the  character  of  the  photorefractive  index  change  occurring  in  our  experiment 
the  electro-optic  coefficient,  T4\,  is  first  determined  by  measuring  the  fringe  shift  with  application 
of  an  applied  field  without  the  focused  self-trapped  beam  present.  The  fringes  shift  one  quarter 
of  a  fringe  when  the  applied  field  Eo  is  switched  from  5.4kV/cm  to  -5.4kV/cm.  This 
measurement  gives  a  refractive  index  change  equal  to  2.5  10'^  and  corresponds  to  an  electro¬ 
optic  coefficient,  r4i  =  -1.42  10''^pmW,  which  is  in  accordance  with  the  literature.'^  In  fig.l  the 
fringe  contours,  with  and  without  the  focused  self-trapped  beam  present  are  shown  we  can  see 
that  the  focusing  or  waveguiding  effect  is  induced  by  an  increase  of  the  refractive  index  in  the 
center  of  the  beam  while  a  strong  decrease  (lO’"')  of  the  index  appears  on  one  side  of  the  beam. 
This  amazingly  strong  decrease  of  the  index  implies  that  a  large  photorefractive  space  charge 
field  is  present.  We  calculate  this  space  charge  field  to  be  about  50kV/cm  when  the  applied  field 
is  only  Eq  =  5kV/cm.  Despite  the  small  value  of  the  electro-optic  coefficient  r4i,  it  is  clearly 
possible  to  induce  large  index  changes.  Physically,  the  large  index  change  is  due  to  the  large 
space  charge  field  that  can  be  created  because  of  both  the  intensity  temperature  resonance  and 
the  low  value  of  the  dielectric  constant  associated  with  InP. 
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THEORETICAL  RESULTS 

In  our  work  we  developed  the  first  model  describing  the  resonant  self-focusing  effects  in 
photorefractive  semiconductors.  The  theory  is  based  on  a  model  that  relies  on  two-donor  levels, 
one  deep  trap  level  and  conduction  of  both  electrons  and  holes.  That  model,  although  successful 
in  explaining  two-wave-mixing,  cannot  explain  self-focusing  effects,  because  it  relies  on  a  single 
periodic  interference  grating  at  low  visibility  (modulation  depth).  Even  a  superposition  of  such 
gratings  cannot  be  applied  to  explain  self-focusing  of  localized  beams.  Yet,  the  underlying 
assumptions  about  the  rate  equations  hold  for  a  localized  beam  and  serve  as  a  starting  point  for 
our  theory. 

We  start  with  the  standard  set  of  equations  in  temporal  steady  state:  the  continuity 
equations  for  electrons  and  holes  respectively,  the  rate  equation  for  deep  traps  level,  the  transport 
equation  and  Gauss  law. 


-  rn^P,  +-T^=  0 
dx 

(la) 

dJ 

^pPt~YpPni+  ^  -  0 

(16) 

+  Y„np,  -t-  -YpP^t  =  0 

(Ic) 

Jn  =  + 

ax 

{\d) 

Jp  =  fJPPpE-qD  ^ 
dx 

(Ic) 

dE  e 

-^-[Nd-Na^p-n-nA 
dx  € 

(1/) 

where  e„  =ef  +^^^„I(x)  and  Cp  =e''’  +(ZipI(x).  ef  and  e'^;  are  the  thermal  emission  rates. 
and  (z5p  are  the  photoionization  cross  sections,  and  are  the  photoionization  cross 
sections,  I  is  the  optical  intensity  and  Nd,  Na  are  the  donors  and  acceptors  concentrations,  nt  and 
Pt  are  the  concentration  of  filled  and  ionized  traps  respectively.  We  seek  an  expression  for  the 
space  charge  field  as  a  function  of  the  optical  intensity,  /,  for  a  single  optical  beam  whose  width 
is  much  larger  than  the  optical  wavelength.  In  such  cases,  the  diffusion  currents  are  negligible. 
Furthermore,  for  typical  value  of  the  parameters,  the  optical  excitation  of  electrons,  as  well  as  the 
thermal  excitation  of  holes  can  be  neglected.  It  is  possible  to  show  (analytically)  that,  within  a 
very  good  approximation,  the  electron  concentration  is  proportional  to  the  holes  concentration. 
The  physical  meaning  of  this  finding  is  that  the  electron  and  hole  currents  are  both  constants,  not 
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only  their  sum.  Keeping  in  mind  that  diffusion  currents  are  neglected,  one  can  easily  see  why 
p(x)=constant  n(x)  (where  the  constant  is  a  positive  number  calculated  at  the  boundary  of  the 
crystal).  Using  this  relation  in  Eq.(lc)  and  substituting  (Id)  in  (Ic),  we  obtain 


E(x)  =  Eo--^^ 

Ires(x)-I(x) 


(2) 


^hexQlns(x)  =  e„n,(x)/((f>p  pt(x))  is  the  local  resonance  intensity.  In  the  same  manner  that 
is  the  rate  of  trap  filling,  (^pP,I,,,  is  the  rate  of  thermal  ionization. 

Therefore  Ir,s(x)  -l(x)  is  proportional  to  the  net  ionization  rate  of  the  trap  level.  In  steady  state, 
this  rate  must  be  compensated  by  the  net  recombination  rate,  which  is  given  by  the  numerator  of 
(2)  times  the  electrons  (or  holes)  concentration.  Note  that,  through  Gauss  law,  1^^^  is  a  function 
of  the  derivative  of  the  electric  field  (the  free  carriers  density  in  Gauss  law  can  be  neglected).  If 
dE/dx  is  small  enough,  is  roughly  constant,  and  thus  the  space  charge  field  can  be  written  as 


E(x)  =  £.  I;  ,  (3) 

I res()  ^  l(x) 

where  Ecm)  =  er,nio/{<l>p  pto)  and  n,Q=Nd-Na.  Clearly,  the  space  charge  field  exhibits 
resonance  effects  when  I-  This  expression  is  no  longer  valid  at  the  resonance,  yet 
numerical  solutions  of  (2)  show  that  it  is  valid  up  to  =0.8/^g^O’  where  l„ax  is  the  peak 


Fig  3:  Soliton  wavefunctions  (a)  and  existence  curves  (b)  for  kerr  (Ancc  I),  saturable  (  A«  oc  /  /(I  +  /)  ) 
and  InP  (  Aw  ocl/(l  -  /)  )  nonlinearities. 
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intensity  of  the  beam.  From  E(x)  of  Eq.  (3)  one  can  find  the  nonlinear  change  in  the  refractive 
index,  An(x)  oc  E(x) ,  substitute  it  in  the  nonlinear  wave  equation  describing  beam 
propagation,^’'*  and  seek  stationary  solutions,  that  is,  spatial  solitons.  Using  techniques  similar  to 
those  of  Refs.  3,4,  we  obtain  the  solitons  wavefunctions  and  the  soliton  existence  curve  shown  in 
Fig.  1. 

The  wavefunction  of  a  soliton  at  -  O.S/^g^o  is  shown  in  Fig  la.  For  comparison,  we 
also  show  the  respective  soliton  wavefunctions  for  the  Kerr  and  the  saturable  nonlinearities.  In 
Fig  lb  we  show  the  existence  curve  of  each  nonlinearity.  As  shown  Fig.  la,  the  normalized 
width  of  the  resonant  (InP)  soliton  is  considerably  narrower  than  that  of  the  Kerr  and  the 
saturable  nonlinearities,  implying  that  solitons  of  the  same  width  require  a  lower  change  in 
An(x)  for  the  resonant  solitons  (a  factor  of  2  at/^^^  =^-^^reso)-  This  is  why  photorefractive 
solitons  in  InP  can  be  observed  in  spite  of  the  fact  that  the  electrooptic  coefficient  of  InP  is  more 
than  100  times  smaller  than  that  of  SEN.  Furthermore,  the  resonant  nonlinearity  of  solitons  in 
InP  does  saturate  (as  the  nonlinearity  of  all  screening  solitons  does),  but  rather  deforms.  As 
the  light  intensity  increases,  the  waveguide  becomes  more  and  more  asymmetric.  This 
asymmetry  stems  from  the  electric  field  derivative  in  (If),  which  is  no  longer  negligible.  In  Fig 
2  we  can  see  a  numerical  solution  for  =  1 .3/,.^.,,, . 


Fig  4:  The  space  charge  field  E(x),  (solid  curve)  and  local  resonant  intensity  Ecsfx)  (dashed)  for  a 
given  optical  intensity  distribution  J(x)  (dash  dotted) .  Note  that  at  the  peak  ,  I  <  Ircs-lrcso  • 

The  resonance  intensity  is  modified  by  the  field  derivative  in  such  a  way  that  the  peak  of  the 
field  is  being  shifted  from  the  center.  It  is  interesting  to  note  that  at  extremum  points  (E’=0) 
expression  (2)  still  holds.  This  leads  to  important  conclusion:  Since  the  numerator  of  (2)  is 
positive  and  the  field  maximum  is  also  positive,  then  the  numerator  must  be  positive  as  well.  In 
other  words  the  maximum  value  of  the  electric  field  will  take  place  at  a  point  for  which 
/(x)  This  explains  why  for  the  beam  undergoes  linear  diffraction':  the 

induced  waveguide  and  the  beam  no  longer  overlap.  These  results  compare  well  with  Fig.  1 . 


6 


While  our  theory  has  been  strong  we  have  also  carried  out  experiments  that  demonstrate 
that  photorefractive  spatial  solitons  can  be  formed  in  semi-insulating  InP:Fe'’^  through  the 
formation  of  self-induced  waveguides.  Our  most  recent  work,  however,  carried  out  under 
support  of  this  grant,  is  on  the  coherent  collision  between  two  solitons  in  InP  at  the 
telecommunications  wavelength  of  1 .3  pm  to  form  a  y-junction. 

For  the  experiment,  light  from  a  diode  laser  is  split  into  two  beams.  Lenses  are  used  to 
independently  focus  each  of  the  beams  onto  the  crystal’s  entrance  face  when  observing  the 
collision  between  two-dimensional  solitons,  while  spherical  lenses  are  used  to  observe  the 
collision  between  two-dimensional  solitons.  The  light  is  polarized  vertically,  along  the  <11 0> 
crystal  axis,  while  the  electric  field  is  applied  horizontally  along  the  <00 1>  axis.  The  beams  are 
arranged  so  that  they  are  parallel  in  the  horizontal  plane.  The  relative  phase  between  the  two 
beams  is  controlled  through  the  use  of  a  piezoelectric  mirror  in  the  one  beam. 

Before  observing  a  collision  in 
InP,  each  beam  is  checked  independently 
to  verify  that  it  forms  a  soliton.  Both 
beams  are  then  launched  together  and 
the  collision  between  one-dimensional 
solitons  is  observed  (Fig  5).  Fig.  2a  Exit  Face 

Ms>i  iiial  DiJl I ticiton 

shows  the  input  of  the  individual  beams  46!imwirfth 

while  Fig.  2b  shows  the  normally 
diffracted  output  of  each  beam.  With  an 
applied  field  of  negative  11.5  kV,  the 
output  diameter  of  each  separately  is 
reduced  (Fig.  2c).  The  collision  of  the 
two  solitons,  when  the  beams  are  in 
phase  (Fig.  2d),  results  in  a  single  peak 
with  a  width  of  23  pm  that  is  located 
centrally  between  the  locations  of  the 
individual  solitons.  We  found  similar 

Fig.  5  Collision  between  two  in  phase  solitons  in  InP 

results  with  a  180-degree  phase 

difference  but  in  this  case  the  two  beam  defected  each  other.  Finally  we  also  investigate  and 
demonstrated  the  formation  of  a  dark  soliton  in  InP.  We  used  the  dark  notch  formed  to  guide  a 
second  laser  beam  (Fig.  6). 

In  summary,  we  have  explained  how  a  mechanism  based  on  conduction  of  both 
(thermal)  electrons  and  (photo-excited)  holes  leads  to  a  resonant  enhancement  (by  more  than  an 
order  of  magnitude)  of  the  space  charge  field.  This  enhancement  is  crucial  to  the  observation  of 


7 


solitons  in  semiconductors  photorefractive  semiconductors  described  below.  An  analytical 
expression  for  the  field  in  the  sub-resonance  regime  was  derived,  and  the  characteristics  of 
solitons  (wavefunctions,  existence  curve)  were  calculated.  Other  observed  effects  such  as  soliton 
collisons  could  be  explained  on  the  same  basis.  Experimentally,  we  have  investigated  the 
collision  between  solitons  and  found  that  the  interaction  between  solitons  is  dependent  upon  the 
relative  phase  between  the  two  beams.  For  a  collision  between  two  solitons  with  the  beams  in- 
phase,  the  solitons  merge  to  form  a  y-junction.  When  the  two  beams  are  out  of  phase,  the  two 
solitons  deflect  each  other.  We  also  demonstated  a  dark  soliton  and  its  ability  to  act  as  a 
waveguide. 
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